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Abstract. A Large Ion Collider Experiment (ALICE) at the Large Hadron Collider (LHC) 
has been built in order to study the Quark-Gluon Plasma (QGP) created in high-energy nu- 
clear collisions. As heavy-flavor quarks are produced at the early stage of the collision, they 
serve as sensitive probes for the QGP. The ALICE detector with its capabilities such as particle 
identification, secondary vertexing and tracking in a high multiplicity environment can address, 
among other measurements, the heavy-flavor sector in heavy-ion collisions. 
We present latest results on the measurement of the nuclear modification factor of open heavy- 
flavors as well as on the measurement of open heavy-flavor azimuthal anisotropy «2 in Pb-Pb 
collisions at ^/s NN =2.76 TeV. 

Open charmed hadrons are reconstructed in the hadronic decay channels D° — > K~7r + , 
D + — > K~7T + 7r + , and D* + — > D°7r + applying a secondary decay- vertex topology. Complemen- 
tary measurements are performed by detecting electrons (muons) from semi-leptonic decays of 
open heavy-flavor hadrons in the central (forward) rapidity region. 



1. Introduction 

A Large Ion Collider Experiment (ALICE) [1] at the Large Hadron Collider (LHC) has been 
built to investigate the properties of a color-deconfined state of strongly interacting matter, that 
is formed in ultra-relativistic Pb-Pb collisions. Heavy-flavor quarks (charm and beauty) are 
well-suited probes to study the hot and dense matter, called Quark-Gluon Plasma (QGP). Due 
to their large mass heavy-flavor quarks are produced in initial hard scattering processes on a 
very short time scale of the order of r ~ l/2mq. After production they interact with the formed 
QGP and because of their early production they are sensitive to the full history of the collision. 
The experimental measurements to access the heavy-flavor sector, which are covered by these 
proceedings, are: 

• Fully reconstructed open charmed mesons D°, D + and D* + in the following hadronic decay 
channels: 

• D° — > K _ 7r + and the charge conjugates. 

• D + —7- K~7r + 7r + and the charge conjugates. 

• D* + — > D°7r + and the charge conjugates. 

• Electrons and muons from charm and beauty semi-leptonic hadron decays. 

One important and well established quantity to characterize the interaction of heavy quarks with 
the medium is the nuclear modification factor Raa studied as a function of transverse momentum 



pr- This quantity is defined as the ratio of the pt spectrum measured in A A collisions to the 
proton-proton spectrum at the same center-of-mass energy yfs scaled by the number N co \\ of 
binary nucleon-nucleon collisions in the AA collision: 

„ / \_ 1 dN AA /dp T _ 1 dN AA /dp T , . 

AA [PT j - (iVcou) • dNpp/dpT ~ {Taa) ■ dapp/dpT > 

with (Taa) being the average nuclear overlap function calculated in a Glauber model of the AA 
collision geometry. N co n and {T AA ) are dependent on the impact parameter b of the colliding 
nuclei, i.e. on the collision centrality. The ratio R aa (pt) equals unity in case of no nuclear 
effects. In a QGP quarks and gluons are expected to interact with the hot medium and lose 
energy. Due to their smaller color charge than gluons, quarks are expected to lose less energy 
than gluons. Furthermore, QCD-theory predicts that the "dead-cone-effect" reduces the energy 
loss of massive charm and beauty quarks with respect to light quarks, see [2H3]- Therefore the 
expected nuclear modification ordering in the range px ^ 10 GeV/c, where the heavy-quark 
masses are not negligible with respect to their momenta, is the following: R AA < R^ A < -^AA' 
with R AA denoting the nuclear modification factor of mostly gluon-originated pions, R^ A the 
nuclear modification factor of open charm D mesons and R^ A the nuclear modification factor 
of open beauty B mesons, see [5]. 

Another important quantity is the azimuthal anisotropy parameter v 2 - In heavy-ion collisions 
the spatial anisotropy in the overlapping zone in non-central collisions is converted into a 
momentum anisotropy of the produced particles, if the in-medium mean free path allows 
for sufficient re-scattering. In this case, the azimuthal distribution of the emitted particles 
reflects the initial anisotropy, which can be parametrized via the coefficients of a Fourier 
expansion 

~ = ^ {1 + 2ui cos(0 - ¥ r ) + 2v 2 cos(20 - 2^ r ) + ...} (2) 
d(f> 2ir 

where *$> r is the azimuth of the reaction plane in each expansion order n. The azimuth of the 
reaction plane for each harmonic n is determined experimentally event-by-event. The Fourier 
coefficients are given by vi,V2,--,v n . The origin of a non-zero v 2 is expected to come from 
different pressure gradients in-reaction plane and out-of-reaction plane at low pr> which is 
known as elliptic flow, and from the difference in the path length which results in different 
energy loss in-reaction plane and out-of-reaction plane at high px- The ^-extraction method 
presented in short here is the " event plane method" . More details about the azimuthal anisotropy 
parameter v 2 , alternative methods of its determination and physics background can be found 
in [6]. 

The proceedings are organized as follows. In Sec. [2] the ALICE detector and the data used 
for the heavy-flavor analysis are described. Section [3] presents the results of the heavy-flavor 
measurements in proton-proton collisions. Next, the i?AA results for heavy-flavor are presented 
in Sec. [4} Finally, Sec. [5] summarizes the ALICE measurement on the azimuthal anisotropy 
parameter v 2 of heavy-flavor hadrons. 



2. The ALICE detector 

For the results presented in the next section the signals of the following sub-detectors of ALICE 
were used: 

• Inner Tracking System: 

The Inner Tracking System (ITS) is located at midrapidity and at radii between 4 cm and 
44 cm around the interaction point and consists of silicon pixel, silicon drift and silicon 
strip detectors. It provides a high spatial resolution, secondary vertexing capabilities and 
is used for trigger purposes. 



• Time Projection Chamber: 

The Time Projection Chamber (TPC) is located at midrapidity at radii between 85 cm 
and 250 cm around the interaction point and is the main tracking device of the ALICE 
detector. 

• Transition Radiation Detector: 

The Transition Radiation Detector (TRD) covers the radial range from 290 cm to 370 cm 
and is located at midrapidity. So far it is used for electron identification at high px only in 
pp collisions for the analyses presented here. 

• Time Of Flight: 

The Time Of Flight system is placed at a radial distance of approximately 400 cm at 
midrapidity and measures the time of flight of particles from the interaction point. The 
time information is used for particle identification. 

• Muon Spectrometer: 

The Muon Spetrometer is placed at forward rapidity (—4.0 < rj < —2.5). A front absorber 
with a thickness of ten interaction lengths A/ and an additional muon filter of 7 A; are used 
to separate muons partially from the background. The remaining background is subtracted 
offline during data analysis, which is described in Sec. |3| 

• VZERO: 

The VZERO detector is composed of two arrays of scintillator tiles covering the full azimuth 
in the pseudo-rapidity regions 2.8 < rj < 5.1 (VZERO- A) and -3.7 < ry < -1.7 (VZERO- 
C). It is used for trigger purposes and collision centrality determination. 

For the results, that are presented in the following, experimental data from the years 2010 and 
2011 were used. A minimum-bias trigger using signals in the two scintillator tiles of the VZERO- 
detector and the silicon pixel detector of the ITS was applied. Moreover, for the Pb-Pb data 
taken in November 2011, a collision centrality trigger was used in order to enhance the number 
of events in semi-central (10-50% most central collisions) and central (0-7.5% most central 
collisions) collisions. This trigger was based on the sum of amplitudes in the VZERO-detector. 
In proton-proton collisions about 100-180 million (depending on the analysis) minimum-bias 
triggered events were analyzed. In Pb-Pb collisions 17 million minimum-bias triggered events 
from the 2010 run and 20 million centrality-triggered events (10-50% most central collisions) 
from the 2011 run were analyzed. 

3. Heavy-flavor production in pp collisions at yfs = 7 TeV and yfs = 2.76 TeV 
In order to determine the heavy-flavor nuclear modification factor a cross-section measurement 
in proton-proton collisions is needed, see Eq. [T] In this section the analysis details of the heavy- 
flavor cross-section in proton-proton collisions are presented. 

Muons from semi-leptonic decays of B mesons and D mesons are measured at forward rapidities 
with ALICE (—4.0 < rj < —2.5) at a/s = 2.76 TeV. The raw muon spectrum is corrected for 
detector efficiency and acceptance using full Monte-Carlo simulations. Moreover, there are two 
main background sources, which contaminate the heavy-flavor muon spectrum even after online 
background suppression (for details see [7j), and need to be subtracted: 

• Muons stemming from primary pion and kaon decays. 

• Muons that are produced inside the absorber due to interaction of the primary kaons and 
pions with the material. 

Monte-Carlo event generators (PHOJET [8j and PYTHIA [91-11 lj) are used to simulate these 
background sources. The kaon and pion spectra generated according to these generators have 
been checked against proton-proton data measured by ALICE. For more details see [7J. The 
background subtracted px-differential heavy-flavor muon spectrum du/dpr taken from [7J is 




Figure 1 . px-differential production cross sections of heavy-flavor decay muons taken from [7] in 
—4.0 < 7} < —2.5 at \fs = 2.76 TeV (left) compared to FONLL calculations [12]. pp-differential 
heavy-flavor electron cross section taken from [13] in \y\ < 0.5 at ^/s = 7 TeV (right, upper 
plot) and the ratio data to FONLL [13 [H] (right, lower plot). 



shown in Fig. [I] (left). The corresponding pQCD prediction [12] agrees with our data within 
uncertainties. 

At midrapidity the heavy-flavor sector is accessed via semi-electronic decays of B mesons and 
D mesons. The crucial ingredient for this measurement is a clean electron sample. For this 
purpose the signals of the TPC, TOF and TRD are combined. After exploiting the signals from 
the three detectors, the remaining hadron contamination is of the order of 10%. More details 
about the particle identification procedure are given in [13J. The main background source for the 
heavy-flavor electron spectrum are electrons originating from photon conversions in the detector 
material and in the beam pipe. The origin of these photons are neutral pion decays, with a 
branching ratio B.R. (ir° — > 77) = 98.8%, abundantly produced in proton-proton collisions at 
midrapidity. A large fraction of this background source is suppressed by requiring a signal from 
the electron-candidate track in the innermost layer of the ITS, which is located at a radial 
distance of only 4 cm from the beam pipe. In order to get the heavy-flavor pr-differential 
electron spectrum, the raw spectrum after all the selections described above is corrected for 
efficiency and acceptance using Monte-Carlo simulations with full detector description. Finally, 
one subtracts the "cocktail" of electrons from background sources from the efficiency corrected 
inclusive electron spectrum. The most contributing components of the cocktail are 

• conversions of non-direct (e.g. from neutral pions) and direct pQCD photons in the beam 
pipe and the innermost silicon pixel layer of the ITS. 



• Dalitz decays of 7r , rj, p, uj and cj) mesons. 

For more details about the cocktail background calculation see [13]. Figure [T] (right, upper 
plot) taken from [13J shows the pr-differential heavy-flavor electron cross section in \y\ < 0.5 at 

= 7 TeV and the ratio data to FONLL [HI 03] (right, lower plot). 
The third method to access the heavy-flavor sector is the full reconstruction of charm mesons 
at midrapidity in the hadronic decay channels. The decays D° — > K~7r + , D + — > K~tt + tt + and 
D* + — > D°7r + are identified via their reconstructed invariant mass. A typical signature of a 
charm meson decay is the presence of two tracks with opposite charges and with a minimum 
distance of approach to the primary vertex not compatible with zero. Once the secondary vertex 
is reconstructed, cut variables can be calculated to reduce the combinatorial background. One 
of the most powerful cuts is the cut on the pointing angle, which is the angle between the 
reconstructed D meson momentum vector and the vector connecting the primary and secondary 
vertex. Another important cut parameter is the decay length, i.e. the distance between the 
primary and secondary vertices. For D mesons this distance is in the order of several 100 //m. 
Further background rejection is achieved by identifying pions and kaons using combined signals, 
which the D meson daughters induce in the TPC and the TOF detectors. The raw yields are 
corrected for efficiency using Monte-Carlo simulation with full detector description. Moreover, 
the efficiency-corrected D meson spectra are corrected for B meson feed-down. This B-feed- 
down correction is needed because according to FONLL (fixed-order next-to-leading log) |12j 
about 15% of the D mesons originate from B meson decays. The B meson feed-down correction 
is performed using the beauty production cross section predicted by FONLL calculation |12j 
and full detector simulation. More details about the D meson measurements in proton-proton 
collisions with ALICE can be found in |15j . The D°, D + and D* + px-differential production 
cross sections taken from [15] in \y\ < 0.5 are shown in Fig. [2j Theoretical predictions based on 
pQCD calculations, namely FONLL O El [16] and GM-VFNS pH HH] are in agreement with 
the data. 

In order to use the measured spectra of fully reconstructed D mesons or electrons originating from 
semi-leptonic D meson and B meson decays at yfs = 7 TeV as references for the corresponding 




Figure 2. pT-differential production cross sections of D° (left), D + (center) and D* + (right) 
at ^fs = 7 TeV in \y\ < 0.5 taken from |15j . compared to pQCD calculations, namely 
FONLL [121 O [E] and GM-VFNS pUCEH]. 

Raa measurement, they need to be scaled down from ^/s = 7 TeV to yfs = 2.76 TeV. The 
scaling factors were obtained as the ratio of the corresponding spectra provided from the pQCD 
FONLL calculations at y/s = 2.76 TeV divided by those at y/s = 7 TeV. For more details 



about the scaling procedure see [19 . The pT-differential cross sections of fully reconstructed 
D mesons were also obtained at the reference center-of-mass energy of \fs = 2.76 TeV, but 
are with limited statistics. It was found that this measurement is consistent with the scaled 
reference at yfs = 7 TeV. For details see [20J. 

4. The heavy-flavor nuclear modification factor Raa 

The selection and background subtraction procedures for the single muon and electron heavy- 
flavor spectrum are the same for Pb-Pb collision data except for minor details described below. 
In case of the muon heavy-flavor spectrum the Pb-Pb analysis differs from the proton-proton 
analysis in treatment of the background subtraction. The spectra of the charged pions and kaons 
used for background subtraction are taken from simulation (see Sec. [3]), while in the Pb-Pb 
collisions the ones measured by ALICE at midrapidity are utilized. Next, they are extrapolated 
to forward rapidities using the rapidity shape from event generators and -Raa information based 
on ATLAS measurements [7] . The Raa for kaons and pions is assumed to be equal at midrapidity 
and forward rapidity but is varied within ± 100% to evaluate a systematic uncertainty for the 
Raa of heavy-flavor decay muons. The nuclear modification factor for the heavy-flavor decay 
muons taken from |7J is presented in Fig. [3] (left) for two different centrality classes, namely 
0-10% and 40-80% most central collisions. A suppression factor of about 3 is observed for 0- 
10% most central collisions and is significantly smaller in more peripheral (40-80% most central 
collisions) collisions. No px-dependence of the -Raa is observed within uncertainties. 
In case of the electron heavy-flavor spectra the difference in the Pb-Pb analysis with respect to 
the proton-proton analysis is that the main input to the cocktail ingredient is not the measured 
neutral pion spectrum but the charged pion spectrum. Moreover, the analysis is performed 
without the TRD and the measurement is restricted to pr < 6 GeV/c to keep the hadron 
contamination below 15%. The nuclear modification factor for the cocktail-subtracted electrons 
is presented in Fig. [3] (right) for two different centrality classes, namely 0-10% and 60-80% most 
central collisions. The large systematic uncertainty is dominated by the electron identification 
uncertainty. A suppression of a factor 1.5-4 is observed for pr > 4 GeV/c in the 0-10% most 
central collisions. In more peripheral collisions (60-80% most central collisions) the suppression 
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Figure 3. (left) -Raa for heavy-flavor decay muons at forward rapidity in central (0-10% most 
central collisions) and peripheral (40-80% most central collisions) Pb-Pb collisions taken from [7j. 
(right) Nuclear modification factors for background-subtracted electrons at central rapidity and 
two different centrality classes, namely 0-10% and 60-80% most central. Statistical (bars) and 
systematic (boxes) uncertainties are shown. 



factor tends to go to unity. 



In case of the fully reconstructed open charm D meson analysis in Pb-Pb collisions also the same 
selection strategy is applied as for the proton-proton case, see Sec. |3j The difference between 
the two analyses is that the topological cuts are chosen in general tighter to deal with the larger 
combinatorial background in the Pb-Pb environment. Moreover, in the Pb-Pb case, a new 
systematic error source is added to the B meson feed-down correction which is not present in 
the proton-proton case. This systematic source comes from the fact that the ratio -Raa/-^aa' 
which is needed for the B meson feed-down correction, is unknown. To estimate the systematic 
uncertainty it is varied in the range 0.3 < -Raa/-^aa < 3. This assumed range represents the 
spread of models predictions for the charm and beauty -Raa [H E] and results on the -Raa for 
J/vt'-mesons stemming from B meson decays measured by the CMS collaboration [21]. This 
variation yields a systematic uncertainty in the final -R AA °^ U P ^° 30%> at high pp- For more 
details about the analysis in Pb-Pb, see [22]. The nuclear modification factor of prompt D°, D + 
and D* + mesons in central (0-20% most central collisions) and peripheral (40-80% most central 
collisions) Pb-Pb collisions taken from [22] is shown in Fig. [4j A strong suppression up to a 
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Figure 4. Nuclear modification factors of D°, D + and D* + mesons for two different centrality 
classes, namely 0-20% (left) and 40-80% (right) most central collisions taken from |22j . 
Statistical (bars) and systematic (empty boxes) uncertainties are shown. 



factor of 3-4 for > 5 GeV/c for all three meson species is observed in central collisions. At 
lower and especially in peripheral collisions the suppression gets significantly weaker. Figure[5] 
(left) taken from [22j shows a comparison of the -Raa averaged over the three meson species and 
a theoretical parametrization of the D meson Raa, which was calculated using a pQCD next-to- 
leading order (NLO) |23j approach with CTEQ6M parton distribution functions [24] including 
only nuclear modified PDF parametrized by EPS09NLO [25J but no heavy-quark energy loss. 
For pt > 6 GeV/c the effect on the Raa caused by shadowing is only ±15% indicating that 
the strong suppression of D mesons is a final-state effect. The analysis of the p-Pb collisions 
data planned to be taken at the beginning of 2013 will enable us to disentangle initial-state 
and final-state effects. Figure [5] (right) taken from [22] shows a comparison between the Raa 
averaged over the three meson species, the -Raa of charged particles [26] and the -Raa of Ja- 
mesons from B meson decays measured by the CMS collaboration |21] . There is an indication for 
-Rcharged < -R^a- However, the indication is not significant enough with the current uncertainties. 
The B meson -Raa measured by the CMS collaboration is clearly larger than the -R A h ? rged but 
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Figure 5. Average nuclear modification factor of D mesons in central Pb-Pb collisions (0— 
20% most central collisions) compared to next-to-leading order (NLO) |23j calculations with 
CTEQ6M parton distribution functions |24j including nuclear modified PDF parametrized by 
EPS09NLO |25j (left) as well as compared to nuclear modification factors of charged particles [26J 
and of J/^-mesons stemming from B meson decays measured by the CMS-collaboration |21j 
(right). Both figures have been taken from [22] . 



due to the large error bars and the coarse binning in no conclusion can be drawn for the 
moment from the comparison of with -R^A- 

5. The heavy-flavor azimuthal anisotropy parameter V2 

In this section the ^"differential azimuthal anisotropy parameter V2 is presented for fully 
reconstructed D° and D + mesons. In order to extract the ^"differential D meson V2, the 
"event plane method" |6] is applied. Due to limited statistics the D meson yield is extracted 
only in two angular regions with respect to the estimated event plane for every event and in 
several pp-bins. The azimuth of the event plane is determined event by event exploiting the 
azimuthal anisotropy of tracks (see [6]) measured by the TPC only. One can show using Eq. [2] 
that the observed D meson V2 is given by 

/ \ _ R obs / n _ R 71 Nin-plane (pr) ~ N out -of-plane (Pt) ,g\ 

* ' ' " ' 4N in 

-plane {Pt) + N out _ of _ plane (p T ) 

with Ni n _ p i ane being the D meson yield extracted in the event plane and N out _of-piane the D meson 
yield extracted perpendicular to the event plane. Moreover, a correction factor R2 accounting 
for the resolution of the event plane [6] is applied. Figure [6] (left) shows the px-differential D°- 
V2 and D + -V2 compared with the charged particle V2 measured by ALICE in the same rapidity 
region [27] . Within the large uncertainties the D meson V2 is consistent with the V2 of charged 
particles. Figure [6] (right) shows the comparison of D meson V2 with calculations using two 
different theoretical models, BAMPS |28] and Aichelin et al. [29]. The model BAMPS [28] is 
a partonic transport model computed with a Boltzmann approach to multi-parton scattering. 
The full space-time evolution of the QGP is obtained by solving the Boltzmann equation for 
on-shell partons and pQCD interactions. The model introduced in |29| is based on a collisional 
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Figure 6. The azimuthal anisotropy parameter V2 of D mesons compared with the one of 
charged particles measured by ALICE [27| (left) and theoretical model calculations, BAMPS [28 j 
and Aichelin et al. [29] (right). Statistical (bars), systematic (empty boxes) uncertainties and 
systematic uncertainties from B-feed-down correction (shaded boxes) are shown. 



energy loss mechanism for heavy quarks. Both models describe the pp-dependence of the D 
meson V2 properly within the large experimental uncertainties. 

6. Conclusion 

ALICE has measured the heavy-flavor nuclear modification factor .Raa in Pb-Pb collisions via 
full reconstruction of open charm mesons and via semi-leptonic decays of charm and beauty 
hadrons. Both in the reconstruction analysis of open charm mesons and in the semi-leptonic 
decay channel of open heavy-flavor mesons a strong suppression comparable to the suppression 
of charged particles up to a factor of 4 is seen in the in central Pb-Pb collisions. A hint of 
^.charged ^ i?^ A is observed in the data. The Raa analysis of the high-luminosity run from 
November 2011 will enable ALICE to precisely measure the heavy-flavor quark -Raa- A first 
analysis of the 2011 run focusing on the pp-differential azimuthal anisotropy parameter V2 reveals 
a preliminary result of the V2 of open charm hadrons. Within the large uncertainties the V2 results 
are properly described by two theoretical models. Finally, a planned p-Pb run at the beginning 
of 2013 will enable ALICE to quantify initial-state nuclear effects. 
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